Oscillating elastic defects: competition and frustration 
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We consider a dynamical generalization of the Eshelby problem: the strain profile due to an inclu- 
sion or "defect" in an isotropic elastic medium. We show that the higher the oscillation frequency 
of the defect, the more localized is the strain field around the defect. We then demonstrate that 
the qualitative nature of the interaction between two defects is strongly dependent on separation, 
frequency and direction, changing from "ferromagnetic" to "antiferromagnetic" like behavior. We 
generalize to a finite density of defects and show that the interactions in assemblies of defects can 
be mapped to XY spin-like models, and describe implications for frustration and frequency-driven 
pattern transitions. 

PACS numbers: 62.20.Dc, 61.72.Qq 



In his seminal 1957 paper Eshelby derived the 
strain fields created by an inhomogencous ellipsoidal 
inclusion in an isotropic elastic medium. This result is 
a cornerstone of the theory of inhomogencous elastic 
media, now routinely used in the physical and engi- 
neering sciences. The result is a statement of how a 
distortion is accommodated in the host material and 
implies that a local perturbation induces long range 
strain fields, slowly decaying as l/r d in d dimensions 
(d > 2). Eshclby's work considers static inhomogeneities 
or strain "defects" only. However, for many applications 
in physics and materials science, we are interested in 
the cooperative behavior of inhomogeneities in which 
the strain is varying or oscillating in time with a 
given frequency We find that this situation is, as a 
function of the defect density and oscillation frequency, 
inherently frustrated, resulting in self-organization of 
the patterns, competing ground states, and sensitivity 
to internal and external perturbations. Such "dynamic" 
defects arise as small polarons in directionally-bonded 
transition metal oxides, including high-temperature 
superconductors, colossal magnetoresistance materials 
and ferroelectrics 0- The collective behavior of these 
polarons in a crystal undergoing distortions, with their 
coupling to charge, spin or polarization, is believed 
to determine the overall macroscopic response. Our 
work also has ramifications for the non-destructive 
evaluation of elastic media. Methods in this field are 
typically based on the vibrational response from a 
defect-free crystal. Here we characterize the behavior of 
oscillating defects (external oscillatory fields inducing 
specific defect patterns and responses) extending the 
conventional analysis to describe the response of elastic 
media in the presence of such defects. 

The dynamical generalization of the Eshelby problem 
has been investigated in the context of engineering sci- 
ences for spherical inclusions Q , and very recently for in- 
clusions of various shapes Q. However, such studies have 
focused on evaluating displacement fields as solutions to 
numerical boundary value problems. Our objective here 
is to understand the effect of dynamics on the nature of 



the elastic interaction itself and its influence on the collec- 
tive behavior of assemblies of defects, using the formalism 
developed in Ref. Q , which allows for tractable analytic 
calculations. We consider the dynamical Eshelby prob- 
lem for localized oscillating defects in two dimensions for 
simplicity, and show that, although the strain fields still 
decay as 1/r 2 far from the defect, the frequency funda- 
mentally affects the nature of deformation. As expected, 
the higher the frequency, the more localized is the defor- 
mation. This renders the interaction between two defects 
strongly frequency (and direction) dependent, but the 
very nature of the interaction changes from "ferromag- 
netic" to "antiferromagnetic" like behavior as a function 
of separation and frequency. We subsequently general- 
ize our results to a finite density of defects. This allows 
us to demonstrate the implications for frequency-driven 
patterning transitions and phase locking in assemblies 
of defects by mapping the elastic interaction energy be- 
tween defects into XY spin-like models with competing 
interactions. 

We use a strain only representation The state of 
strain is defined by three fields e* related to the displace- 
ments along the x-axis (u) and the y-axis (v) as follows: 
ei = (u x + v y )/VZ, e 2 = (u y +v x )/\/2, e 3 = (u x — v y )/V2 
(the subscripts x and y indicate differentiation). We 
write an elastic energy which includes gradient terms: 
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where A\ and A are the bulk and shear moduli of the 
isotropic elastic material, respectively, and g is a strain 
gradient coefficient. Assuming first an overdamped dy- 
namics (we comment later on the undcrdamped case), 
the equations to be solved read, in real space: 
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-(A + gV 2 )e 2 + 2dxdyX , (3) 

-(A + g\7 2 )e 3 + (d 2 x-d 2 y)X , (4) 

V 2 ei — 2dxdye 2 + {d 2 x — d 2 y)e 3 . (5) 
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The last equation is the compatibility condition, reflect- 
ing the fact that the three strain fields are not inde- 
pendent; A is a Lagrange multiplier enforcing this con- 
straint. The boundary conditions are taken to be peri- 
odic in space. We now add n oscillating defects localized 
around the positions (fi)i—i... n . We assume for speci- 
ficity that the oscillation is in e 3 . Considering strictly 
point defects, that is 63(7^, t) = eoS(r — fl) sm(ujot + ipi), 
with wo and ipi denoting oscillation frequency and phase, 
would lead to unphysical logarithmic divergences of the 
Green functions for the system of equations 10)- ©• 
Thus, we enforce the regularization: J ei(fl,t)g a (r) dr = 
J e 2 {fi,t)g a {r) dr = 0, / e 3 (fl,t)g a (r) dr = e sin(u> i + 
ifi), where g a {r) = exp(— r 2 /2cr)/(27rer), and a is cho- 
sen small enough that the defect is physically localized. 
We have checked that our results, which apply to the far 
field created by the defects, do not qualitatively depend 
on the regularization. 

We start with the case of one defect. As the problem 
is linear, we can construct the solution as a superposition 
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of (fr± , the elementary solutions of the problem, defined 
as follows: 



£3) 



2 v 



-(A x 
-(A + gW 
-(A + gW 2 



jl(H) 

'</4 

4 12) - 

4 13) - 



- V 2 A + . 9(T (r)e 
2dxdyX , 
(d 2 x - d 2 y)X , 



±tuJot , (6) 
(7) 
(8) 



The 



and 



are solutions of the same set of equa- 



tions with the oscillatory excitation in the second and 
third equations, respectively. The 4>± can be analyti- 
cally calculated in Fourier space. Writing the ej as linear 
combinations of the 4>± , and finding the coefficients by 
enforcing the regularized oscillating defect conditions, we 
obtain the expressions for the fields. We give here explic- 
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Two important results can be deduced from this exprcs- 
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sion. First, due to the compatibility condition, <fy + is 

not continuous around k = 0. This creates the same 
1/r 2 tails as in the static case. The Eshelby result thus 
extends to oscillating defects. Second, from Eq. I|1U|) we 
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see that the larger u>q, the more spread out is 4>\ (k). 
In real space this implies that the larger luq, the more 
localized is the deformation created by the defect. 
This will be of primary importance for the interactions 
between defects, as detailed below. 

These results are illustrated in Figs.n an d|2 showing 
strain profiles for different uiq. A comment on the un- 
derdamped case is in order: an undcrdamped dynamics 
would not remove the discontinuity at k = created by 




FIG. 1: Profile of the strain field e% along (Ox) for frequencies 
Wo = 1 (blue diamonds), 10 (red circles), 100 (black crosses). 
The smaller u>o, the wider the profile. The parameters used 
are A\ = 4, A — 3, g — 3, a — 0.01. Inset shows the surface 
plot of e3. Notice the anisotropy of the field. 
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FIG. 2: Log-Log plot of the strain field e% along the diagonal 
for too = 1, 10, 100 (same symbols as in Fig.0. We have also 
added 1/r 2 fits as guides to the eye. The parameters used are 



At = 4, A = 3, g 



0.01. 



the compatibility equation; thus, the 1/r 2 decay is also 
valid in this case. The qualitative effect of increasing loq 
would not be modified either, although there would be 
some quantitative differences from the overdamped case. 
Finally, we have focused here on defects created by a 
locally oscillating e 3 strain; the solutions for locally os- 
cillating e\, e2 strains, or combinations of the three strain 
components, can also be obtained, and are qualitatively 
similar. There is one exception to this statement: the 

functions (f>^ (k) are continuous around k = and even 
infinitely diffcrentiable. Their tail in real space is thus 
exponential instead of power law. This implies that the 
e\ strain field created by a local e\ defect in an isotropic 
clastic medium decays exponentially away from the de- 
fect with rate po; if the defect is oscillating, the expo- 

1 /2 

ncntial decay rate p(u>o) grows with ujq as iv . To our 
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FIG. 3: U(ip) for two defects along the x-axis. The distance 
between the two defects is d — 2 (blue curves), d — 4 (red 
curves) or d = 8 (black curves). The frequency is u>o = 1 
(solid lines), 10 (dashed lines) or 100 (filled circles). The 
parameters used are Ai — 4, A = 3, g — 3, a = 0.01. 




FIG. 4: U(<p) for two defects along the diagonal. The distance 
between the two defects is d = 1.96 (blue curves), d = 3.92 
(red curves) or d = 7.84 (black curves). The frequency is 
wo = 1 (solid lines), 10 (dashed lines) or 100 (filled circles). 
The parameters used are A\ = 4, A = 3, g — 3, a — 0.01. 



knowledge, this particular case has not been emphasized 
in the literature; it would be interesting to realize its 
experimental signatures. 

We now consider two defects oscillating with the same 
frequency u>q, with the goal of studying the interactions 
between them. Let the two defects be centered at f = 0, 
with phase ip\ = 0, and r — fo, with phase ip2 = p. We 
express the fields as linear combinations of <f>± (r) and 

<j)^\r — fo)- As above, the twelve complex coefficients 
are found by enforcing the regularized oscillating defects 
conditions. For each strain configuration, it is easy to 
calculate the elastic energy stored in the system from 
Eq. Q. Averaging this energy over one oscillation pe- 
riod, we obtain the interaction energy of the two defects 
U((p,fo,u)o). Similar calculations in the static case were 
performed, for instance, by Eshelby in Ref. 0. 

We consider first that the two defects are pinned, and 
study the function U(tp) at fixed fo and loq. It turns out 
that the energy can be approximated by an XY spin inter- 
action term, U (ip,f o,uo) = — J(rb, <jJq) cos ip. A positive 
J then corresponds to a "ferromagnetic", phase-locking 
interaction, and a negative J to an "antiferromagnetic" 
one. As expected, J decreases to zero at large distances; 
from the previous one-defect calculations, it can be an- 
ticipated that at fixed distance, J also approaches zero 
in the large u>o limit. These effects are seen in Fig. [3] 
which shows the energy as a function of tp for two defects 
along the x-axis. This figure also emphasizes a striking 
effect: J may also change sign with varying fo or luq. For 
small enough too and small enough distance, the interac- 
tion is ferromagnetic, and the defects tend to phase-lock; 
for larger distances and luq, the interaction is antiferro- 
magnetic, and the energy is minimized for a maximum 
phase difference tp = n. Since the fields are anisotropic, 
the picture is different for defects situated along the diag- 
onal; in this case, the interaction is always ferromagnetic, 
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FIG. 5: The solid line is the approximate boundary be- 
tween "antiferromagnetic" and more complex ground states 
for chains of defects parallel to the x or y axis (as predicted by 
our two-defect calculation). The dashed line is the crossover 
line between "ferromagnetic" and "antiferromagnetic" near- 
est neighbor interactions. The insets show the results of a 
numerical energy minimization for three different cases. The 
parameters used are A\ = 4, A — 3, g = 3, a = 0.01. 

see Fig. 0] 

This qualitative change in the interactions between 
defects from varying the frequency or the distance has 
important consequences when considering the collective 
properties of assemblies of defects, as shown below. 
We have focused here on oscillating strain fields 
for specificity; for oscillating e± or e-i strain fields, the 
results are quantitatively different, but the main feature 
remains: changes in the frequency or the distance have 
dramatic effects on defect interactions. 

Using a superposition of the elementary solutions 
(jr± \ it is possible to construct the strain profiles asso- 



4 



dated with an assembly of N defects. The regularized 
oscillating conditions have now to be enforced for each 
defect, which yields the 6N coefficients of the linear 
combination. Calculating the energy (averaged over one 
period) as a function of all the phases, we are now in a 
position to study the collective behavior of the assembly 
of defects. For the sake of specificity, we focus again on 
oscillating e% strains. As a working example, we consider 
chains of equally spaced defects, along the x-axis or the 
diagonal, with different spacings and frequencies, and 
study the ground state of the system as a function of 
the phases of all defects. 

To obtain a qualitative understanding of the pattern- 
ing, it is useful to infer the behavior of N defects from 
the pairwise interactions studied above. Note however 
that summing pairwise interactions is an approximation, 
neglecting three-body effects and beyond. These effects 
are included in the full numerical calculations. When 
defects are aligned along a diagonal, the pairwise 
interactions are always phase-locking, see Fig. there 
is no frustration and the ground state should be a 
perfectly phase-locked, or "ferromagnetic" state, for 
all lattice spacings and frequencies. This is indeed 
observed in our full calculations. The case of defects 
aligned along the x-axis is richer, see Fig. 01 As "anti- 
ferromagnetic" interactions come into play, frustration 
and competition between different types of interactions, 
such as long-period patterns found in ANNNI spin 
systems Q, are likely to develop. However, at large 
enough spacing along the chain and/or large enough 
frequency, the interaction between nearest neighbors is 
dominant, and "antiferromagnetic" ; the ground state 
is then expected to show neighboring defects with a 
phase difference of n. At small lattice spacing and 
frequency, the interaction between nearest neighbors is 
"ferromagnetic" ; the interaction between next nearest 
neighbors however is "antiferromagnetic" . Frustration 
will then develop, and a non-trivial ground state may be 
expected. At intermediate lattice spacing and frequency, 
all interactions are "antiferromagnetic" , but the nearest 
neighbor interactions do not dominate; in this case 
also, a complicated ground state should be expected. 
Although this picture based on two-body interactions is 
approximate, it provides qualitatively correct results, see 
Fig. [S] Figure [S] also shows the numerically determined 
boundary between "antiferromagnetic" ground states 
(where nearest neighbor interactions dominate) and 
more complicated ones, as expected from the two-defect 



calculations. The boundary (dashed line) between 
ferromagnetic and antiferromagnetic nearest neighbor 
interactions is also sketched. We have performed our 
calculations here for special arrangements of defects, 
but these results demonstrate the existence of collective 
behavior, controlled by the frequency or interparticle 
distance. 

In conclusion, we have generalized the well known 
Eshelby study of static elastic inhomogeneities Q, |(| to 
dynamical oscillating defects, focusing on the qualitative 
properties of the resulting strain fields. Our main results 
for a single defect are: a f/r 2 decay of the strain fields 
far from the defect, and that the higher the frequency, 
the more localized is the strain perturbation. These 
results have important consequences for the interaction 
between oscillating defects, which slowly decays at large 
distances, and is suppressed by increasing ujq. A more 
detailed study of the two oscillating defects situation 
shows a more dramatic effect: varying the frequency, 
or the distance, can result in a qualitative change in 
the interaction, from "ferromagnetic" (phase-locking) 
to "antiferromagnetic" (anti-phasc-locking). Our work 
raises the possibility of controlling collective patterning 
of defects by the frequency or the defect density; we 
have explicitly demonstrated such effects by performing 
N-defect calculations. 

Applications of our work include non-destructive 
evaluation of clastic media with internal oscillating 
defects, and the collective behavior of small-polaron 
dopant sites in dircctionally bonded electronic materi- 
als @- There are natural extensions of our results to 
the dynamics of multiple defects. Thus, phase- locking 
among defects can be studied if specific phase dynamics 
controlling the relaxation pathways is added. Again, 
if the center-of-mass of the defects is given dynamics, 
we can expect self-assembly of defect patterns driven 
by the induced elastic fields; since the interaction range 
is reduced for higher frequency defect oscillations, this 
will reduce the domain of defect ordering. When there 
are competing ground states due to ferromagnetic and 
antiferromagnetic interactions, we also expect multiscale 
"glassy" dynamics. These extensions of our results will 
be reported elsewhere. 
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